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INTRODUCTION 


This report presents the initial study on a guidance system for improved 
accuracy projectiles that employs three, pre-emplaced beacons, as shown 
in figure 1. 


Three beacons, B,, Bz, and B;, are emplaced by some means in the 
vicinity of a fixed target, T. The locations of the target and the beacons 
are assumed known ‘a priori', within specified statistical errors, from 
intelligence sources. A transponder at each beacon site transmits signals 
to a receiver in the projectile, where a micro-processor calculates the 
ranges, R,, Re, and R;, from the projectile to the beacons, B,, Bz, and 
B,, respectively. 


Utilizing the a priori beacon, target location data and the measured 
ranges, the on-board micro-processor calculates the instantaneous pro- 
jectile position and the projectile-to-target vector. This vector is then 
utilized as the steering signal input to a navigator which generates tre 
control forces necessary to guide the projectile to the target. 


The objectives of this study are: 


1. Determine the feasibility of the three-beacon guidance system 
to locate, precisely, the coordinates of the projectile relative to an inertial 
reference system, given the beacon and target locations and the projectile- 
to-beacon ranges. 


2. Determine the needed beacon and target location accuracies 
needed to achieve a specified circular error probability (CEP) at impact. 


3. Determine the feasibility and applicability of using the projected 
planar guidance system, as utilized by Sandia Corporation (ref. 1), for 
solving this problem. 


GEOMETRY OF THE THREE-BEACON SYSTEM 


The three beacons, B,, Bz, and Bs, are assumed to be emplaced at 
some nominally known locations with respect to an inertial reference 
frame (x, y, Zz) as shown in figure 2. 


A beacon-based coordinate system (x°, y’, z’) will be erected with the 
origin, O’, at the location of beacon B,, the x°-axis along the line from B, 
to B,, all three beacons in the x’-y” plane, and the z°-axis orthogonal to 
the x’°-y’ plane. 


In order to assign identities to the pre-emplaced beacons so that the 
geometry shown in figure 2 is accurate, the following methodology is 
employed: 


1. Beacon B, is chosen as that beacon closest in distance to the 
origin of the inertia! reference frame. 


2. A tentative choice for B, is made from the remaining two beacons 
and the vector cross product B, B, x B, B; is determined and then normalized 


to determine the unit vector, k’, in the z’ direction. 


3. This process is repeated after interchanging the identities of 
beacons B, and B,. 


4. The scalar product k-k’ is then calculated for each k’. 


5. Beacon B, is then taken as that choice of B, that results ina 
positive scalar product, k:k’, with B,; as the remaining beacon. 


This identity assignment will be accomplished at the firing site, or at some 
ground station, once the beacon locations have been determined. 


With respect to the inertial reference frame, the assumed coordinates of 
the beacons are: 


Ros, =j a, +j b, tke, (1) 
ao g a a 

‘op, ag +j be +kce 

Sop, = ag +j bs + kc; 


The equation of a plane passing through al! three beacons is: 
(x-a,) (y-b,) (Z-c,) (2) 
(a,-a,) (b,-b,) (c,-c,)} =0 
|(a3~a,) (bs ~b, ) (cs -C,) 
Upon expanding equation (2), we have 
A (x-a,) +B (y-b,) +C (z-c,) =90 (3) 


or 
Ax +By +Cz=D (4) 


where 
A = (b,-b,) (c3-c,) - (bg-b,) (c,-c,) (5) 
B = (ag-a,) (ce-c,) - (ag-a,) (Cg -c,) 
C = (a,-a,) (b3-b,) - (a3-a,) (b,-b,) 
D =Aa, + Bb, + Cc, 


We can normalize equation (4) to obtain the equation of the plane 
through B,, Bz, Bs as, 


xX COS G, + y COS Bg + Z COS Y3 =P (6) 
where 
cos ad, =A/(A? +B? +2)? (7) 
2 


cos B, = B/ (A? + B* +C?) 
cos ¥, =C/(A? + B? +C?) 
p =D/(A? +B? +C?)? 
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The three cosines are the direction cosines of the plane and p is the 
perpendicular distance from the origin of the inertial reference frame to 
the plane. Thus, a unit vector normal to this plane is 


k’ = (cos a3) i + (cos B,) j + (cos y3) k (8) 
The unit vector i’ along B, B, is obtained as follows: 
BBs = OBy OB, ell day ats) +7 (Or BY) BR ee.) (9) 
Normalizing this we obtain 
i’ = (cosa,) i + (cos B,) j + (cos y;) k (10) 
where 
cos a, = (ag: --a,)/ [lag --a,)* + (bg - by)* + ( - c,)?12 (11) 
cos B, = (b, - b,)/[ (ag - a)? + (be, - b,)? + (cy -¢,)?]}? 


cos V1 = (co = c,)/[ (a, a + (be an b,)? + (co ~¢,)2]? 


The unit vector fie can now be obtained from the cross product 


j° =k’ xi’ (12) 
” al 3 j k 
J) =] cos ag, cos Bs, COS ¥3 (13) 
COS Qy cos Bp, cos ¥j 


which yields 


j° = (cos ag) i + (cos Bo) i 4 (cos Y,) k (14) 
where 

COS @, = cos B; cos y,; ~ cos B, cos yy (15) 

cos B, = COS a, COS Y3 - COS @, cos y, 


cos = COS ds, COS ~— COS @, COS 
2 3 1 1 3 


The relationship between the inertial and beacon coordinate systems 
can be expressed as: 


cosa, cosfB, cosy, (16) 


coOS @ cos PB, COSY, 


| a a | 


ins 
j’ 
ik’ 


cOS G3 cosB3; cos ¥3 


This matrix is orthogonal, so that the inverse matrix is equal to the 
transpose. Therefore, 


COS Q, COSQ, COS @z (17) 


i i 
j cos B, cosP, cos Bs; i‘ 
c COSY; COS ¥_ COS ¥3 k’ 


It is now desirable to express the coordinates of the beacons in the 
beacon coordinate system (x°, y', z’) as functions of the beacon coordinates 
in the inertial system. From our selected geometry we have 

B, (a,, b,, C1); By (a2, De, Co); Bs (ay, Ds, C3) (18) 
in the inertial coordinate frame, while 

B, (0, 0, 0); Bz (ag, 0, 0,); Bs (as, bs, 0) (19) 
in the beacon coordinate system. Thus, 


<2 —> —_ _7 = 2 = 
B,B, = OB, - OB, = dg i ° = i (a, ~ a,) + j (b, - b,) +k (c, > €,) (20) 


taking the scalar product of this and i’ obtained from equation (16) we have 
the result, 


(ag i’) 57 =a, = (a, -a,)? + (bz - b,)* + (cy = G5) : (21) 


—s Z » 
Similarly, taking the scalar product of B,B; with both i’ andj” from 
equation (16) we have 


a ae p>.» AG > Par, , 
(B,B;) -i° = (a, i +bgj ) +i =as (22) 
—2 7? A>. rae 2 zy , 
(B,B,;) -j° = (a3 i +bg3 j ) - Jj = bs 
or 
a3 = (a, a a,) cos Qy + (bs b, ) cos B, + (cy a C,) cos Y1 (23) 
b, = (as - a,) cos a, + (bg - b,) cos Bg + (cg ~ c,) cos ¥2 


COORDINATES OF THE PROJECTILE 


At any instant of time during the flight, let the coordinates of the pro- 
jectile be (x, y, z) in the inertial coordinate system, and (x°, y’, Z) in 
the beacon coordinate system. Similarly, the target locations at any time are 
(ay, by, cy) and (ap bt c +), as shown in figure 3. The time of arrival 
(TOA) or distance measuring equipment (DME) measured distances from the 
projectile to beacons B,, Bz, and Bs are R,, Rz, and R;, respectively. 


These ranges can be expressed in terms of the projectile position as: 


R,? = (x’)? + (y’)? + (z’)? (24) 
Ry? = (x’ - ag)? + (y’)? + (27)? 
Rs” = (x’ - ag)? + (y’ - by)? + (z’)? 


The projectile can be located by a sequential! solution of the foregoing 
equations in the following manner: 


R,* = Re = 2as x” = fe, )* (25) 
or, 
x’ =| Ry? -R,? + @3)* | / ag (26) 
Similarly, 
R,? -R,? = 2a, x ° - (af)? +2 bg y’ - (bs)? (27) 
or 
y= [Rt =Ry" (a, )° (b= 2a, x” | 126; (28) 


and finally, 
2° = 4VfRy? - (x)? - (y’)? (29) 


With the projectile location now known in the beacon coordinate system, 
we note that 


—» ——»> —— > 
P =OB, +B,P (30) 
or | 
xit+tyj+zk =(a, i +b, j te, Ky 2! fey’ Jo ea Ke) (31) 


Substituting equation (16) into equation (31) and solving for the inertial 
coordinates of the projectile, we have 
x=a, +x’ cosa, ty’ COS ad, +z COS a, (32) 
y=b, +x’ cos B; + y’ cos By +z’ cos By 
Z=C, +x’ cosy, +y COS ¥, +z COSY; 
the inverse of which is 
x” = (x ~a,) cosa, + (y — bj) cos B, * (z- ¢,) cosy, (33) 
y’ = (x-a,) cos ag + (y — by) cos B, + (z —cq) cos y2 
7” 


= (x-a,) cos ag + (y - b,) cos B, + (z-c,) cos ¥3 


Now, the projectile to target vector PT can be shown to be 


— >) > z " Fe: 
PT = OP - OT = (x - a7) i + (y - by) j + (2-1) k (34) 
or 
PT =O’P - O'T = (x’ - at) i’ + (y’ - bp) i + (2 - ef) K SP) 


The results of either equation (34) or (35) can be utilized for guidance 
purposes. In this report, equation (34) is used in order to simplify the 
projectile equations of motion. 


It should be noted at this time, that equation (29) is ambiguous with 
respect to the projectile altitude in the beacon coordinate system. This 
result is readily understandable when we consider that the solution to the 
projectile location problem consists of finding the loci or intersection of 
three non-colinear spheres with different radii, which, for the general case, 
results in two discrete points. 


To resolve this ambiguity, consider when the projectile is at time t = to, 
the start of active guidance. At this time, the altitude in the inertial system 
was assumed to be roughly the apogee of the ballistic trajectory. At this 
point, we can estimate the projectile position (x, y, z) in the inertial system 
from a Simple ballistic solution, and since we know, a priori, the values of 
a,, by, C,;, COS a3, cos B3, and cos 73, we can determine, with the use of 
equation (33), the algebraic sign of the projectile altitude in the beacon 
coordinate system. in addition, by continuity arguments, we can state that 
if z’ decreases and ultimately changes sign (fig. 4) it does so smoothly, so 
that the sign of the calculated z° would change subsequent to every zero 
crossing. 


TRAJECTORY CONSIDERATIONS AND PROJECTED PLANAR GUIDANCE 
SCHEME 


It has been shown that the three-beacon range measuring system has 
the capability of determining the projectile coordinates and the projectile- 
to-target vector at any instant of time, in either coordinate system. There- 
fore, given perfect range measurements and true beacon and target locations, 
the achievable CEP at impact depends on the chosen guidance law coupled 
with the projectile dynamics and the altitudes at the start and termination of 
guidance. 


The simulation developed in subsequent sections permits a comparative 
evaluation of different navigation schemes for a point-mass projectile, with 
provisions for modeling bias-type errors in beacon and target locations and 
random errors in range measurements. 


Sandia (ref. 1) has reported that a "projected planar guidance" scheme 
is superior to a proportional navigation system when used with pre-emplaced 
beacons. This guidance scheme can be explained with reference to figure 5. 


The projectile is constrained to be steered in the x — y plane only, 
utilizing the command signals (x(t) - ay), and (y(t) - br). Thus, the 
velocity of the projectile in the x - y plane becomes 


Vp=S- [x(t) - ap) i +S fy(t - byl j (36) 
or 
mae = ssa 
Vp = x(t) i + y(t)j (37) 


In the projected planar guidance scheme, the planar projectile velocity 
is smoothly decreased during the guidance phase so that at the termination 
of guidance, t=t,, | 


x(te) =O (38) 


x (tg) =o 
y(te) =o 
y (te) =0 


X (te) - ay =O 
y (te) - by =0 


At termination of guidance, the projectile is at a point directly above the 
target, with zero planar velocity and acceleration. For the remainder of the 
trajectory, the projectile falls vertically with respect to the x - y inertial 
plane, until target impact. 


Equation (38) indicates that this scheme makes no use of the altitude 
information, which can be provided by the three-beacon system. So, in 
essence, a data item is being discarded by the scheme despite the chance 
that its inclusion may result in improved guidance accuracy. 


THE PROJECTILE EQUATIONS OF MOTION AND THE PROJECTED 
PLANAR GUIDANCE LAW 


Projectile Equations of Motion 
A flat-earth, three-degrees-of-freedom point mass model was assumed 
for the projectile with drag, Fp, gravity, mg, and applied control forces, 


oe Fy, and F,, as shown in figure 6. 


From Newton's Law, 


F=ma (39) 
where 

F=f that Ret thy) Fz k (40) 
and 

=) d? = d? 

a=ae “ae Few; (41) 
or 

a=xX ity j+tzk (42) 


Thus, the gravity force ts, 


Fg =~ (mg) k (43) 
and the drag force is 


Fp = - 4 pV? Cp A uy 
where 


U ee ou Flap - 2h J wk (44) 
Vv Vv V Vv 
so that 


Fp=-ip CpAV [xityj+zk| (45) 


Substituting equations (43) and (45) into equation (40) and (42) into 
equation (39) and solving for the component accelerations, we have, 


on ep . (46) 
x =— -kp— Vx 
m m 
F CHA 
ee 
par inp My 
: CyA 
hae gp Vz 


G92 + (y)? + | : 
a 
22000 


< 
T 
a 


where m = Projectile mass 


A = Projectile cross-sectional area 
g = Acceleration of gravity 
p = Air density at altitude z 
Pas= Air density at sea level 
Cp = Coefficient of drag 
er Fy, rz = Applied control forces 
‘ x, Y, Z = Projectile component velocities 


x, Y, Z = Projectile component accelerations 


The Navigation Law 


The particular guidance scheme studied in this report was the projected 
planar guidance scheme as reported by Sandia (Ref. 1). In this scheme, 
F, is assumed to be zero, while Fx and Fy, are each linearly proportional] 
to a combination of their respective position and velocity commands. This 
results in 


x ; 
we Kp Gy # Ky &y (47) 
F 
ee ns Y 
vn Kp ey + Ky ey 
F =0 
Z 
Where 
ey = x(t) - a7 (48) 
ey = y(t) - by 


The equations for the three-beacon system, as previously developed, 
provide e, and e,, directly, but do not provide their derivatives. The 
required derivatives can be approximated by a lead-lag network that has 
a transfer function of the form 


KYS (49) 
S 


-— 
sal 


where A is selected to be large, when compared with the frequency content 
of the command signals, e, and ey - 


Combining equation (49) and equation (47) we obtain the transfer 
function of the navigation controller 


F K,,S 
=a Kp + . &x ae! 
eS 5 
F K.,S 
ar wy, e 
m= \Sp* y 
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SIMULATION PROCEDURE 


The method of computer simulation is shown in block diagram form 
in figure 7. 


Initial Computations 


At the start of the simulation, the true beacon and true target 
coordinates, 


C; PS1ewz2, By "Tt 


with respect to the inertial reference frame are specified, as well ~s the 
bias error in beacon and target locations, 


(aj), =aj + Aa. 

(bj) ¢ = b; + Abj i=1, 2, By F 

(cp)¢ = cj + Ac; 
The bias errors, Aaj, Ab;j, Acj, are modeled as independent, Gaussian 
random variables, with zero means and specified variances. These bias 
errors are obtained at the start of each simulation and are held constant 


for the duration of the simulation. 


Using both the true and biased values, the true and calculated direction 
cosines are determined from equations (7), (11), and (15) 


cos Gj, cos B;, cos Yj m= 1, 2,3 
(cos aj)c, (cos Bj)c, (cosyj). i=1, 2, 3 


and the true and calculated beacon locations in the beacon reference frame 
are determined from equations (21) and (23) 


az, &, b3 
(ab)o, (ah)o, (bie 
Dynamic Computations 


During the dynamic simulation, the projectile equations of motion, 
equation (46), are integrated to determine the true location of the orojectile 
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in the inertial reference frame (x, y, z). The coordinate transformation, 
given by equation (33), is then utilized to obtain the true projectile position 
with respect to the beacon reference frame (x', y’, z'). The true ranges 
are then determined from equation (24). 


Range measurement errors are then introduced by adding computer 
generated noise to each true range, 


Ri (t) =Ri(t) + ep (t) (Sty 2 3 


Rj 


where each noise generator, Ep (t), is independent and preduces Gaussian 


i 
white noise, with zero mean and identical variance. 


The projectile coordinate calculator then determines the calculated 
projectile position with respect to the beacon reference frame, utilizing 
the measured ranges, Ring: and the calculated beacon locations in the 


beacon reference frame, (a3),., (a3), and (b3), in equations (26), (28), 
and (29). 


The coordinate transformation given by equation (32) is then employed 
to arrive at the calculated projectile position in the inertial reference frame, 
which, along with the calculated target location (a. , Bie 7 cy ), is used 

Te Te Cc 


in equation (48) to determine the calculated steering inputs (ey ; e, ) to 

Cc Cc 
the steering control system. The control system then produces forces F,, 
and F.., according to the relationship given in equation (50), to alter the 
projectile flight dynamics. 


It should be noted that the extent of the computations required by the 
projectile's on-board processor are given by blocks 4 through 6 in figure 7. 
These are relatively simple, straight-forward calculations; however, they 
require the on-board storage of the calculated beacon locations in both the 
inertial and beacon reference frames, the calculated target location in the 
‘nertial frame, and the calculated direction cosines. These data must be 
calculated and furnished to the projectile processor prior to firing or during 
the early portion of the projectile flight prior to the start of active guidance. 
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Extension of the Simulation 


The simulation method, as shown in figure 7, was developed to model 
the fixed beacon reference system. 


The simulation as it exists, however, can be readily adapted to include 
beacons, each with its own proper motion. In this case, the beacon positions 
are allowed to vary with time, relative to the inertial reference frame. Thus, 
ai, Bix Gj, Become aj (t).. bptt), et), =, 2, 3, Wherérap(t).; b; (t) and 
c;(t) are known at any time. Beacon errors can then be viewed as true 
positions plus random errors, rather than as true positions plus bias errors, 
as indicated previously. 


Similarly, the target location which was assumed to be fixed and known 
within a mapping bias error, can assume a proper motion by allowing its 
coordinates to vary with time, a7 (t), b+ (t), c_ (t) , With associated bias 
and random white noise errors. 


tn order to include a moving target scenario in the simulation, some 
mechanism must be postulated for measuring the target position as a function 
of time. Speculation as to the nature of this mechanism, however, is beyond 
the scope of this report. 


STATIC SENSITIVITY AND ERROR ANALYS!S 
The determination of the sensitivity of the impact CEP to beacon 
location errors is a prime objective of this study. The degree of sophisti- 
cation of the guidance !faw is immaterial if the inherent inaccuracies in- 


troduced by the three-beacon guidance system cause the projectile to 
exceed the desired CEP at impact. 


The impact CEP sensitivity of the three-beacon system when only errors 
in beacon location are considered. is investigated. This case is referred to 
as the static case. 


Static CEP Sensitivity 


In considering only beacon location errors, we make the following 
assumptions: 


1. Range measurements are erroriess 


2. Target location is known exactly 
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3. The guidance system is ideal and capable of realizing the projected 
planar guidance trajectory exactly. Thus, the calculated trajectory will 
result in projectile impact at the target. 


Although we know the target location exactly, the steering signals 
Q, = X_ (t) - ay (51) 
ey = Ps ee oa 
depend on the calculated projectile position (denoted by subscript "c"), 
which is in error due to the beacon location errors. Therefore, the 
calculated trajectory impacts at the true target location, T, while the 


true trajectory impacts at some other point, T¢, as shown in figure 8. 


At t = te, the projectile is at some prescribed altitude, Z., and the 
projected planar guidance will have achieved the following conditions: 


Xo = ay Yc = by (52) 
“=o ¥.= 6 
X. = 0 fae 


The constraints imposed by equation (52) can be used to find the actual 
projectile position at t = te as follows using equation (33) 


ew (53) 
x. * (ay - ay) COS Gy. + (b+ - b; 0) cos Bie ze C1) cos Yi. 


? 


Ye 


(ay - ay) COS dp # Aba = bi) cos Bo. + = C1) COS Yo. 


? 


Z. = (a= a, 0) cos Oy 42 (loa = by) cos Bs. + Ze = C0) cos Me. 

The true ranges — the projectile to the beacons can then be found by 
R? = (xc)? + (yc)? + (z¢)? (54) 

RE = (xg - az)? + (yg)? + (20)? 


RS = (xe 7 as)? + (yg > bs)? + (zg)? 
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We can now obtain the actual projectile position in the beacon-based 
coordinate system 


Ri - R$ + (a3)? (55) 
2a5 


R? - R3 + (aj)? + (bg)? - 2ay x’ 
4 = 2b; 
2 SRG Ft SS iy ie 


Thus, the actual projectile position in the inertial reference frame is 


xX = a, +x° cosa, +y’ coSa, +z’ cosa, (56) 
y = b, +x’ cos PB, +y’ cos B, +z cosB, 
Z = ¢, +x’ cosy, +y’ cosy, +Z° cosy; 


The coordinates of the actual projectile impact are given by equation 
(56), if the constraints of equation (52) are satisfied for the true trajectory. 
Since the guidance system calculations are based upon the calculated 
trajectory, the constraints of equation (52) will not, for the general case, 
be satisfied, since x and y at t = te are non-zero quantities. The actual 
impact point of the true trajectory can be computed by integrating the 
ballistic equations of motion for t>tf, using the actual position and velocity 
at t = tp as the initial conditions. 


Rather than pursuing this line of computation, we can obtain more 
easily an approximate determination of the true point of impact. At t= te, 
z is relatively small compared to the maximum ordinate of the trajectory, 
so that it can be assumed that x, y, and z are essentially constant from 
t= te to impact. Due to the geometry of the trajectory, at t = tf 


Z>xXx, y (57) 
We can assume that the total velocity is essentially z(tp). Since the 
calculated trajectory closely approximates the true trajectory, we can state 


that 


Zt;) = Zc (tr) , Z(t-) = Z. (tf) (58) 
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Thus, the approximate time to impact, At, can be calculated as 


we 


Then, the impact point is given by 


x) = x(t) + (te) At (60) 


y, =y (te) + y(tp) At 


and the planar impact error at the target is 


p= (<a) * = Gy = by)" (61) 


The projectile position at t = t¢ is obtained by utilizing equations (53) 
to (56), while the velocities x(t), y(t¢), and z(tf) are obtained in the 
following manner: 

Differentiating equation (53) we have 


Xo = 2% (cos yi)¢ (62) 


Yo= Ze (cos y2). 
Zc= 2. (coS ¥3)¢ 
Next, differentiating equation (54) and utilizing the results of 
equation (62), we have 
2R,R, =2 Zc Xe (cos Wl Fy. (COS Ya)q.+ 7 (cos Ya)c| (63) 
2R,R, =2 dics (x.- ag )(cos ¥1)6 + Ve (COs ¥2)¢ + 2 (cos ¥s)¢ | 
Cc 
RyRy = 228 | (xg - a3.) (cos Ye + (¥e- bs.) (cos Y2)¢ + 2% (cos Ys)c| 
Differentiating equation (55) and utilizing the results of equation (63), 


we have 
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x = (az /a2) Zo (COS ¥1)¢ (€4) 


yY =Ze¢ (b3 /bs) (cos Y2)_— + (1/bs;) fas. 95 (az /az) (cos ¥,)-} 


Z =2Z. =Ze (cos Ve). 


Finally, differentiating equation (56) and using the results of 
equation (64) we have, 


de bs 
x =z cos a, (cos y,). + —_! cosa (cos y,),. + ly Mee) 
ru ? ] 1 ? 2 2 ? 
Cc a» Cc bs < b3 
oS 
(a3. ——a ) cos a, (cos ¥,). + cosa, (cos 73), 
as b; 
ae Cc Cc 1 
Y=Z ; cos B,; (cos y1)¢ + |p | cos By (cos y2), +H. 
da» bs bs 
a; 3. 
(a3. = “a Jcos By (cos ¥1), +cos B3 (cos ¥3). 


Thus we can obtain the planar velocities att = tp in terms of the 
calculated vertical velocity, z-., the true and calculated beacon locations 
in the beacon coordinate system, a5, a3, b3, aa a3 b3 and the 


direction cosines between the beacon and inertial reference frames, cos Kine 
cos B;, cos yj, |= 1, 2, 3. 


Error Statistics 

Beacon Errors 

We have assumed that the components of the beacon location error are 
Ax, Ay, and Az. These component errors are modeled as independent 


Gaussian random variables, with zero mean and equal variance, Of: . The 
spherical beacon error is taken as 


Pp = (Ax)? + (Ay)? + (Az)? (66) 
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For the Gaussian model assumed above, the probability density function 
for pp is 


Pr” pur 
fo (oa) nye Pe exp(= = (67) 
B “FB? 3 2 
T on 2P, . 
This distribution is characterized by 
‘i Z i: 
average {pp} = Pp = a2 o, = 1.60 o. (68) 
and 
CEP {p,} = 1.540, (69) 
where CEP {p,,} is derived from Pr by 
< = 
P Pp CEP (p,,) 0.5 (70) 


Target Impact Errors 
At the time of impact, the coordinates of the impact error are 


(Ax) + =x)- a (71) 


- 
(Ay)_ =y, ~ b+ 


and the planar circular error is 
p= Y (Ax)7+ (Ay)> (72) 


At this point, we will assume that (Ax) + and (Ay)+ can be approxi- 
mated by independent Gaussian random variables, with zero mean and 
equal variance. This implies that p+ Is Rayleigh distributed with a 
probability density function of the form 


2 


P+. P : 
fr (p-.) = aa EXP< = 702 (73) 
T By 
This distribution is characterized by 
average {p+} =p = 1.25 O_. (74) 


and 


CEP {p_} = 1.18 0 (75) 


+ 


Utilizing both equations (74) and (75) we obtain 
CEP (Pr 3 = 0.944 P+ (76) 


Dividing this by equation (69) we obtain 


p p 
CEP UT) ~ 9 613) (77) 


p 
CEP 4 B} O., 


which gives the ratio of target-to-beacon CEP's as a function of average 
impact error and beacon standard deviation. 


For the computer simulations, 0,, was specified, and the Ax, Ay, and 
Az of the beacons were obtained by multiplying Op by the output of a 
random Gaussian distribution subroutine with a zero mean and a variance 
of 1. The average impact error, p__, was then calculated, based ona 
large number of static simulations. 


COMPUTER SIMULATION RESULTS 
Static Simulation 
Program and Parameters 


A computer simulation was developed to imp!ement equations (51) 
through (65). The program was written in FORTRAN Extended, Version 4, 
for use ona CDC 6600 computer system and used a proprietary CDC 
mathematic/statistical routine, NRAND, to generate nermally distributed, 
pseudo-random numbers. 


For the simulations considered, the beacon coordinate system is 

simply translated from the inertial reference system (no rotations involved). 
Thus, the locations of the three beacons and the target, relative ta the 
inertial frame, are taken to be: 

By, §& (RX, RYO) 

By : (RX BS, RY <O) 

B; +: (RX, RY # BS, 0) 

~T ¢ TAT, Bie) 


The inputs to the program are: 
ZF- Zs, guidance cut-off altitude (meters) 
VTERM - Projectile terminal velocity (meters/second) 
SIGB - On: beacon location standard deviation (meters) 


BS - Beacon spacing (kilometer) 


RX, RY - X and Y components of beacon #1 relative to inertial 
frame (meters) 


C1, C2, C3, CT - Altitudes of beacons #1, #2, #3 and the target 
in the inertial reference frame (meters) . 


The program outputs were: 
XBAR - Mean impact error (meters) 
STDEV - Impact error standard deviation (meters) 


RTCEP - Circular error probability (CEP) of the impact points 
about the target location (meters) . 


A listing of the computer program statements, as utilized, is given in 
figure 9. The input variables had the following values: 


ZF - 0,500, 1000, 2000 meters 
VTERM - 0, 200, 250, 300 meters/second 
SIGB - 1, 4, 7, 10, 30, 50 meters 


BS - 1, 2, 3, 4, 5, 10, 15, 20, 25, 30, 35, 40, 45, 50, 55, 
60 kilometers 


RX, RY - 15,000 meters 


C1, C2, C3, CT - 0 meters 
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Simulation Results 


Appendix A presents the results of the static simulations. It should 
be noted that the results presented for each beacon spacing represents 
the statistical results obtained from 100 independent Monte-Carlo simula- 
tions for the indicated input parameters. 


The ratio, CEP {p_.}/CEP {Pp }, as a function of beacon spacing and 
for fixed values of beacon leeation variance and projectile terminal 
velocity, is given in figure 10. It can be seen that for relatively close 
beacon spacing, an Increase in the guidance cut-off altitude, ZF, produces 
a corresponding increase in impact error, which is consistent with the 
results of equation (60). It is interesting to note, however, that there 
exists a value of beacon spacing, (BS),,j,, such that CEP {p}/CEP {pp} 
attains a minimum for all ZF and all BS 4(BS),,jn. Figure 11 is an expanded 
plot of a portion of Figure 10, and indicates that this beacon spacing is on 
the order of 10 kilometers. 


Figure 12 presents CEP {p__}/CEP {p,} as a function of guidance cut- 
off altitude and terminal velocity for fixed-beacon variance and beacon 
spacing. It appears that CEP (pr)/CEP {pp} attains a minimum at different 
guidance cut-off altitudes for each of the terminal velocities considered. It 
can be seen, however, that the impact error over a 2000 meter increase in 
guidance cut-off altitude increases by only 11 percent. We can, therefore, 
state that the impact error asa function of guidance cutoff altitude is 
relatively insensitive to variations in projectile terminal velocity. 


The variation of CEP {p--} as a function of CEP {pg} and terminal 
velocity for a fixed beacon spacing and guidance cut-off altitude is 
presented in figure 13. !t can be seen from this that CEP {p._} varies 
linerary with respect to CEP {pp} and is independent of terminal velocity. 


Dynamic Simulation 
Program and Parameters 


Dynamic simulation of a typical projectile employing projected planar 
guidance and three beacons was implemented on a CDC 6600 computer 
utilizing CDC's Continuous System Simulation Language Ill (CSSL 3). This 
is a computer program designed to facilitate the representation and simula- 
tion of continuous dynamic systems. The language provides simple and 
straight forward programming of problems involving differential equations 
as opposed to an equivalent FORTRAN program. 
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Beacon location errors are obtained from a Gaussian random-number 
CSSL3 subroutine at the start of each simulation and then held constant 
for the duration of the run. Each beacon error is independent of the 
other. 

Ranging errors from the projectile to the beacons are obtained from 
the random number subroutine and are inputted at each computation interval 
during the run. These simulated ranging errors are made independently 
beacon-to~-beacon. 

Figure 14 presents a listing of the CSSL3 simulation program. 

The program inputs are as follows: 

VO - Projectile muzzle velocity (meter/second) 
QE - Quadrant elevation (mils) 

AZ - Firing azimuth (mils) 

A - Projectile cross-sectional area (meters) 


XMASS - Projectile mass (kilograms) 


AT, BT, CT - Location of target relative to gun position 
(meters) 


Al, B1, Cl 

A2, B2, C2 Locations of beacons #1, #2, #3, respectively, 
relative to gun position (meters) 

A3, B3, C3 

SIGMA - Standard deviation of beacon location error 


RSIG - Standard deviation of ranging error 


P, Q - Coefficients of the guidance transfer function 
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The program outputs were: 


xX, Y, Z- Projectile location with respect to the gun 
position (meters) 


T - Projectile flight time (seconds) 


For all simulations, the following parameters were fixed at the values 
indicated: 


Muzzle velocity - 555.2 meters/seconds 
Quadrant elevation - 804.9 mils 
Firing azimuth - 800.0 mils 
Beacon variance - 5.0 meters 
Range variance - 25.0 meters 
Beacon spacing - 5 kilometers 
Twenty simulations were run for each of the following four cases: 
Case | - Control errors only 
Case II - Control and beacon bias errors 
Case III - Control and ranging errors 
Case IV ~ Control, beacon bias, and ranging errors 
and for each of the following two target geometries: 


Geometry | - Target within triangle formed by beacon 
locations. 


Geometry II - Target outside triangle formed by beacon 
locations. 
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Simulation Results 


Appendix B presents the results of the dynamic simulations for each 
error case and target location. The target CEP is calculated using equation 
pecs 


Figure 15 shows the CEP at the target for Case | errors for both 
target geometries. The CEP for this case is quite small, on the order of 
0.44 meters, but not actually zero. This error can be reduced by optimi- 
zing the control system parameters. 


Figure 16 presents a plot of the impact points and CEP about the target 
for error Case Il and Geometry |. The calculated CEP for this case is 7.5 
meters. Figure 17 presents the results for error Case Ill, Geometry |. The 
calculated CEP is 11.2 meters. 


Figure 18 presents the results for error Case IV, Geometry |. The 
calculated CEP is 14.5 meters. 


As can be seen from figures 15 to 18, errors due to beacon bias and 
ranging are roughly comparable and an order of magnitude larger than 
control errors for this particular geometry. It can also be seen that the 
total system CEP is not the summation of the CEP's due to the individual 
errors. That is, 


* | (cep! ~ (cep) *| (cep) 


(p 
CEP Case IV Case Ill Case | Case Il 


= ile ) + (o ) 
CEP Case | CEP Case | 


Although Geometry | is a reasonable case to assume, it does not reflect 
the expected deployment of this type of projectile guidance system. The 
expected deployment of the guidance beacons would have all three beacons 
located to one side of the front edge of the battle area (FEBA) while the 
target is located on the other side of the FEBA. This senario was simulated 
by placing the target outside of the triangle formed by the beacon locations 
and it resulted in Geometry II. 
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Figure 15 and figures 19 through 21 present the results for error 
Cases | through IV with Geometry II. The calculated CEP's for each case 
are given below: 


Case | CEP = 0.44 meters 
Case ll CEP = 14.7 meters 
Case Ill CEP = 13.3 meters 


Case IV CEP = 19.6 meters 


It can be readily seen that the increase in system CEP for Geometry II 
is due mainly to the increase in error contributed by beacon bias. This is 
apparent since small errors in direction cosine computation result in 
larger impact errors in Geometry II than in Geometry | due to the increased 
beacon-to-target ranges encountered in Geometry II. 


CONCLUSIONS 


1. The feasibility of accurately locating a projectile in space utilizing 
a three-beacon guidance system has been shown. 


2. The projected planar guidance scheme as utilized by Sandia (ref. 1) 
appears feasible to be used in conjunction with the three-beacon guidance 
system. 


3. Beacon and target location accuracies of + 5 meters in each axis 
are acceptable. It can be expected that beacon location accuracies on the 
order of + 0.1 meter can be achieved by utilizing standard artillery surveying 
techniques (ref. 2). 


4, Effects of beacon spacing and target location, with respect to the 
beacons, appear to be the most significant parameters effecting the 
achieveable CEP about the target. 


RECOMMENDATIONS 


Further effort should be directed toward evaluating conceptual, 
extended-range projectiles utilizing the three-beacon guidance system. 
Included in this effort should be the choice and optimization of a particular 
control scheme, the determination of an optimum control law, the modifica~ 
tion of the projectile equations of motion to include all known aerodynamic 
effects, and the upgrading of the dynamic simulation from a three to six- 
degree-of-freedom program. 
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Figure 5. Projected planar guidance scheme 
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Figure 6. Forces acting on a point-mass projectile 
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Figure 9. Static simulation computer program 
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Figure 9. (continued) 
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Figure 9. (continued) 
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APPENDIX B 


DYNAMIC SIMULATION RESULTS 
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Run 


Table B1. 


Number 


LO CO™ AN & WP — 


ee ee ee ee ee ee 
WO CONMA OS WH MY — © 


Dynamic simulation results, Case 1, Geometry | 


(p cep), = 0.44 meters 


X - Impact Y - Impact 
9899.19 9899.19 
9899.19 9899.19 


Table B2. Dynamic simulation results, Case Il, Geometry | 


( p cep). = 7.5 meters 


Run Number X- Impact Y- Impact 
] 9897.59 9902.09 
2 9890.23 9894.29 
3 9898 .99 9899.09 
4 9871.49 9887 .93 
5 9895.54 9902.22 
6 9884.54 9895.71 
Fi 9893.07 9902 .38 
8 9897.01 9901.94 
9 9898 .85 9889.61 

10 9904.89 9898 .09 
1] 9899.04 9896 .97 
12 9892.30 9903.00 
13 9894.40 9903.20 
14 9903.70 9897.03 
15 - 9896.71 9900.49 
16 9884.15 9896.71 
17 9904 .80 9903.43 
18 9904.25 9895.42 


1g 9891.66 9905.31 


Table B3. Dynamic simulation results, Case II|, Geometry | 


( p cep). = 11.2 meters 


Run Number X - Impact Y - Impact 
] 9912.64 9906.01 
2 9889. 3] 9899 .82 
3 9893.92 9897.53 
4 9904.52 9903.42 
5 9905.54 9909.07 
6 9889 .06 9900.36 
] 9896 .88 9902.26 
8 9920.22 9909.26 
9 9884.17 9885.28 

10 9897 .48 9910.56 
1] 9884.28 9886.06 
lz 9904.56 9899.83 
13 9907.95 9897.85 
14 9900.02 9920.68 
15 9901.84 9904.36 
16 9896.69 9903.16 
17 9893 .87 9880.06 
18 9896 .54 9883.02 


Table B4. Dynamic simulation results, Case IV, Geometry | 


(P cep), = 14.5 meters 


Run Number X - Impact Y - Impact 
] 9911.19 9909.09 
2 9880.56 9894.11 
3 9893.90 9897.61 
4 9890.77 9906.27 
5 9902.04 9912.29 
6 9888 .44 9910.95 
Fi 9890.95 9905.96 
8 9918.22 9912.22 
fe) 9883.93 9875.90 

10 9903.35 9909.65 
1] 9884.27 9884.01 
12 9897.85 9903.82 
13 9903.33 9902.05 
14 9904.64 9918.62 
15 9899 .53 9891.92 
16 9881.80 9900.85 
17 9899.68 9884.5] 
18 9901.75 9879.37 
19 9888 .99 9903.77 
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Table B5. Dynamic simulation results, Case !, Geometry II 


( Pp cep). = 0.44 meters 


Run Number X- Impact Y- Impact 
9899.19 9899.19 


WON DOP WD 


19 9899.19 9899.19 
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Table B6. Dynamic simulation results, Case II, Geometry II 


( P cep), = 14.7 meters 


Run Number X - Impact Y- Impact 
] 9897.75 9901.25 
2 9872.85 9888. 37 
3 9897.13 9900.95 
4 9900. 4/ 9928.21 
5 9878.70 9893.80 
6 9909.37 9928.70 
7 9925.24 9911.73 
8 9899.12 9908.16 
9 9916.19 9883.02 
10 9911.81 9897.35 
1] 9910.31] 9901.49 
12 9902.41 9905.92 
13 9886.09 9893.88 
14 9890.30 9883.37 
15 9898.23 9896.87 
16 9893.92 9892 .00 
17 . 9890.62 9891.01 
18 9904.41 9910.73 
19 9876.92 9884.41 


Table B7. Dynamic simulation results, Case III, Geometry II 


( p cep). = 13.3 meters 


Run Number X - Impact Y - Impact 
] 9917.77 9910.95 
2 9891.85 9900.03 
3 9895.44 9899.69 
4 9904.86 9906.10 
5 9901.91 9904.92 
6 9891.11 9902.10 
7 9899.64 9905.57 
8 9924.2] 9916.19 
9 9880. 36 9881.50 
10 9903.22 9916.23 
1] 9877.79 9878.58 
12 9905.67 9901.39 
3 9905.4] 9897.70 
14 9908.51 9925.64 
15 9904.11 9906.29 
16 9897.86 9901.45 
17 9888 . 16 9876.50 
18 9896.77 9886.50 


18 9902.83 9903.88 
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Table B8. Dynamic simulation results, Case IV, Geometry II 


(p cep) = 19.6 meters 
C 
Run Number X_- Impact Y - Impact 
] 9916.57 9913.06 
Z 9865.80 9889.35 
3 9893.56 9901.67 
4 9906 .35 9935.33 
5 9881.64 9899.64 
6 9901.45 9931.75 
7 9925.94 9918.28 
8 9924 .32 9925.39 
9 9897.49 9865.53 
10 9916.03 9914.48 
1] 9889.00 9881.01 
12 9909.07 9908 .31 
13 9892 .53 9892.54 
14 9899.68 9909.89 
15 9903.30 9904.16 
16 9892.76 9894.43 
17 9879.81 9868 .55 
18 9902.18 9898.17 
19 9880.73 9889.26 
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